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摘要 : 30~300 MHz 的 甚 高 频 (Very High Frequency, VHF) 频段 是 重要 的 射电 天 文 波段 ， 
该 频段 观测 采用 天 线 阵 组 阵 方 式 。 稀 布 阵 列 具 有 空间 分 辩 率 高 、 副 闪电 平 低 以 及 造价 低 等 优 
点 ， 进 一 步 的 天 线 阵 综 合 加 权 可 以 对 天 线 阵 主办 波束 进行 有 效 赋 形 ， 对 最 大 旁 准 副 闪 
( Maximum Side Lobe) 电 平 和 远 区 机 办 ( Far Side Lobe) 电 平 进行 抑制 。 首 先 回 顾 了 射电 天 文 其 
高 频 稀 布 阵列 研究 发 展 和 现状 ， 以 及 将 会 遇 到 的 难点 ， 提 出 了 首先 优化 最 优 稀 布 天 线 阵 元 排 
布 ， 进 一 步 提 出 融合 高 性 能 计算 平台 +FPGA SOPC 的 稀 布 其 高 频 射 电 天 文 阵列 信号 处 理 结构 
体系 ,在 图 形 处 理 器 ( Graphics Processing Unit, GPU) 或 者 云 计 算 平 台 上 完成 对 天 线 阵 各 阵 元 
频 点 加 权 参 数 的 计算 ， 然 后 通过 高 速 总 线 将 计算 参数 下 发 到 前 端的 信号 处 理 板 中 ， 通 过 
FPGA SOPC 完成 对 加 权 参 数 的 配 发 。 进 一 步 分 析 计 算 了 多 波束 情况 下 的 数据 率 ， 可 以 实现 
实时 的 参数 配置 。 本 文成 果 为 下 一 步 大 规模 甚 高 频 天 线 阵 架设 提供 了 技术 依据 。 
关键 词 : 甚 高 频 ; 优化 算法 ; 稀 布 阵列 ; 阵列 综合 ; 低 副 关 
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甚 高 频 (30~300 MHz) 频 段 是 重要 的 天 文 观测 频段 ， 在 空间 天 气 和 太阳 射电 、 太 阳 系 内 行星 
爆发 、 脉 冲 星 和 和 暂 现 源 ”” 、 宇 宙 早 期 再 电离 "等 诸多 领域 有 重要 的 研究 价值 。 为 此 ,， 近 20 年 诸 
多 射电 天 文 观测 新 设备 均 针 对 该 频段 建立 ， 具 有 代表 性 的 有 VLA( Very Large Array), GMRT( Giant 
Meterwave Radio Telescope), SKA (Square Kilometer Array), LOFAR ( LOw Frequency ARray ) 和 LWA 
(Long Wavelength Array) 46575! 。 此 频段 的 特点 是 观测 信号 波长 长 ， 只 能 采用 多 个 单元 振子 天 线 组 阵 
的 方式 进行 观测 ， 通 过 加 权 形 成 高 增益 波束 实现 高 灵敏 度 。 同 时 根据 射电 天 文 的 特点 ， 其 高 频段 需要 
宽带 宽 观测 和 宽 角度 扫描 。 另 一 方面 ， 该 频段 无 线 电 环境 复杂 ， 需 要 接收 机 系统 兼顾 较 高 的 动态 范 
围 ， 在 条 件 允 许 的 情况 下 可 以 通过 空间 “ 零 陷 ”抑制 来 自 某 些 方向 的 干扰 信号 。 

该 频段 射电 天 文 实现 高 空间 高 时 间 分 辩 率 及 高 灵敏 观测 的 基本 思路 是 由 几 百 个 振子 天 线 ( Antenna- 
Beam) 组 成 一 个 站 ， 利 用 波束 合成 技术 ( Station-Beam) 实现 对 射电 源 的 跟踪 ; 然后 再 由 几 百 甚至 上 千 
个 这 样 的 天 线 站 组 成 阵 ， 每 个 站 的 波束 进行 两 两 相关 ， 得 到 高 空间 分 辨 率 的 图 像 。 具 体 的 波束 合成 过 
程 如 图 1。 上 述 两 个 层次 波束 合成 的 区 别 ， 站 内 波束 是 由 一 维 数据 流 组 成 ， 代 表 信 号 来 自 一 个 固定 的 
天 区 ， 而 站 间 的 相关 干涉 是 创建 天 空 的 二 维 快照 图 。 即 站 内 波束 合成 完成 了 输入 信号 的 权重 相 加 ， 而 
站 间 的 则 完成 输入 信号 的 相干 "7 。 
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图 1 其 高 频 天 线 阵 波束 合成 层次 
Fig.1 Three levels of VHF array beamforming pipeline 
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前 述 多 个 天 线 阵列 是 采用 先 构建 子 阵 (Tle) ， 通 过 子 阵 中 模拟 方法 (Tile beamformer) 构建 一 个 初 
步 的 子 波束 ， 再 通过 子 阵 之 间 进 一 步 的 站 级 别 波束 合成 (Station beamformer) 。 这 样 可 以 降低 一 定 的 数 


据 率 ,但 是 在 天 线 阵 指向 灵活 性 方面 有 所 下 降 。 


平方 公里 阵列 中 的 低频 阵 (SKA-low) 提 出 了 全 数字 化 阵列 ， 即 每 个 阵 元 天 线 的 X\Y 两 个 极 化 均 
连接 模拟 接收 机 和 高 速 数字 采集 。 离 散 傅 里 叶 变 换 ( Discrete Fourier Transformation，DFT) 预 处 理 链 
路 的 构架 如 图 20 2 ， 这 样 可 以 做 到 更 高 的 空间 分 辨 率 ， 不 存在 子 阵 列 合 成 后 波束 变 窗 、 视 场 下 降 等 


问题 。 
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图 2 基于 离散 傅 里 叶 变 换 的 全 数字 化 阵列 示意 图 
Fig.2 All digital array based on DFT 
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1 稀 布 阵列 及 其 在 射电 天 文中 的 优势 


采用 传统 的 规则 布 阵 (Regular Lattice) 实现 超 宽 带 阵 列 综合 面临 的 主要 问题 是 其 在 进行 大 范围 天 
区 扫描 时 会 出 现 机 办， 栅 注 的 存在 严重 制约 站 级 别 波束 的 性 能 。 然 而 ， 由 于 超 宽带 阵列 在 进行 大 范围 
天 区 扫描 时 ， 最 小 阵 元 间距 通常 无 法 满足 上 限 频率 方向 图 的 无 栅 流 条 件 

d < A/(l + sin|0,, |), (1) 

其 中 ，d 为 阵 元 间距 ; 9 为 扫描 角度 ; 在 波长 A 一 定 的 条 件 下 ， 随 着 9 的 增 大 无 法 满足 阵列 方向 图 机 
办 限 制 条 件 。 为 了 解决 栅 办 问题 ， 超 宽带 阵列 需要 设计 成 阵 元 不 等 间隔 分 布 的 稀 布 阵列 ， 即 超 宽 带 稀 
布 阵 列 。 稀 布 阵列 可 以 打破 阵列 辐射 能 量 在 空间 的 周期 性 分 布 ， 达 到 消除 栅 办 的 目的 。 此 外 稀 布 阵列 
还 具有 其 他 一 系列 优点 ， 比 如 可 以 增 大 阵列 口径 ， 提 高 空间 分 辨 率 ， 无 需 幅 度 加 权 即 可 实现 低 副 办。 
为 此 ， 目 前 超 宽带 稀 布 阵列 的 综合 技术 已 成 为 当前 阵列 信号 处 理 领域 的 研究 热点 。 

非 均匀 间隔 的 稀 布 天 线 能 够 节省 大 量 成 本 ， 采 用 天 线 阵 元 在 一 定 孔 径 上 稀 布 的 方法 实现 阵列 设 
计 ， 可 以 在 相同 的 天 线 阵 元 数目 条 件 下 ， 得 到 更 高 的 增益 和 空间 分 辨 率 ， 从 而 简化 结构 ， 降 低 造价 ， 
成 为 天 线 领域 的 研究 热点 2 。 目 前 应 用 多 种 算法 可 以 将 均匀 阵 ( 圆 阵 、 直 线 阵 等 ) 进行 稀 布 布 阵 ， 以 
达到 减少 天 线 阵 元 、 抑 制 阵列 天 线 方向 图 的 旁 辩 电 平 、 降 低 成 本 的 目的 | 。 

借助 上 述 天 线 阵 发 展 的 诸多 成 就 ， 终 端 可 以 快速 处 理 大 规模 天 线 组 成 的 超 宽 带 稀 布 阵列 ， 特 别针 
对 目前 其 高 频频 段 射 电 天 文 天 线 阵 ， 这 将 是 一 种 行 之 有 效 的 研究 手段 。 

在 射电 天 文 方面 ， 针 对 目前 其 高 频频 段 天 线 阵 的 问题 和 压力 ， 在 平方 公里 阵列 等 后 续 阵 列 提出 
时 ， 将 该 阵列 定义 为 稀 布 阵列 配置 ' 汪 ” ， 特 别 是 在 低频 频段 。 文 [24] 分 析 了 密集 阵列 ( Dense array ) 
和 稀 布 阵列 在 灵敏 度 ( Sensitive) 、 视 场 (FOV ) 和 巡天 速度 (Survey speed) 方 面 的 参数 比较 ， 如 表 1。 表 
1 中，Aw 为 紧密 阵列 下 最 小 阵 单元 间隔 (m) , A 为 观测 波长 ，N, 为 天 线 数量 ，Ni 为 接收 通道 数量 。 
首先 在 有 效 面积 方面 ， 稀 布 阵 列 的 最 低频 率 比 密集 阵列 提高 20 倍 ; 其 次 ,将 两 者 的 巡天 速度 做 成 比 
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根据 (2) 式 ， 两 者 的 巡天 速度 比较 如 图 3。 由 图 3 可 以 看 出 ,在 500 MHz 以 下 的 甚 高 频频 段 ， 稀 布 阵 
列 比 密集 阵列 具有 更 大 的 优势 。 


表 1 密集 阵列 和 稀 布 阵列 参数 对 比 Survey speed comparison 
T 
1 o 
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图 3 在 巡天 速度 方面 低频 段 稀 布 阵列 占 优 
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n—J i, ate ER AE FERE RESZ CL) SNB d), BERT 0 增 大 ， 超 宽带 阵列 进 
行 大 角度 扫描 (如 + 上 45 ) 时 ， 无 法 满足 阵列 方向 图 顶 浴 限制 条 件 。 在 扫描 角度 较 大 时 ， 从 上 限 频 率 开 
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始 在 方向 图 上 出 现 栅 办 ， 该 栅 办 的 绝对 电 平 与 主 
办 相近 ， 会 对 观测 数据 带 来 灾难 性 的 影响 。 

目前 ， 稀 布 阵 列 特别 是 超 宽 带 稀 布 阵 列 的 研 
究 ， 通 常 选 取 最 高 旁 兴 电 平 ( Maximum Side Lobe 
Level, MSLL) 抑制 度 为 主要 优化 目标 ， 需 要 满足 
超 宽 带 范围 内 的 全 频段 方向 图 无 栅 兴 的 条 件 。 在 
稀 布 布局 方案 确定 后 ， 对 稀 布 阵列 进行 各 阵 元 加 
权 ， 在 波束 宽度 、 副 闪电 平 、 远 区 栅 为 电 平 多 个 
参数 和 维度 上 进行 赋 形 ， 比 传统 的 布 阵 方式 有 较 
大 的 提高 ， 特 别 是 副 办 电 平 能 抑制 在 30 dB 以 
上 。 图 4 是 文 [24] 分 析 仿 真 的 稀 布 阵列 后 进 一 
步 波束 加 权 赋 形 后 的 波束 情况 。 

LOFAR 在 LBA( 10 ~ 80 MHz) 频段 开始 考虑 
稀 布 阵列 的 优势 ， 为 了 获取 稀 布 因子 ，LOFAR 
采用 随机 指数 分 布 布局 二 ， 如 图 5。 

在 后 续 的 天 线 阵 中 ，LWA 的 设计 者 开始 考 
虑 将 天 线 图 综合 优化 方案 融合 到 稀 布 阵列 中 。 优 
化 目标 主要 有 两 种 : (1) 降低 副 瓣 电 平 ; (2) 增 
加 空间 频 域内 的 采样 点 数 029 。 

在 LWA 中 , LWA-1 的 椭圆 形 阵列 配置 方面 
进行 了 优化 ， 以 最 小 化 最 高 副 斩 为 优化 目标 ， 优 
化 函数 为 天 线 阵 接收 矢量 ”3 ， 


He Sah e 
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天 线 方向 图 函数 作为 优化 目标 ， 和 迭代 多 次 ， 
最 终 优 化 阵 元 的 布局 位 置 如 图 6。 优化 目标 仅 为 
天 线 阵 指向 天 顶 时 ， 降 低 最 高 副 办 电 平 ， 采 用 的 
方案 是 将 256 阵 元 初始 化 布局 于 一 个 标准 的 六 边 
形 阵列 中 ， 通 过 逐步 迭代 调整 各 个 阵 元 的 位 置 ， 
并 以 该 天 线 阵 站 内 区 域 面积 和 最 小 可 容忍 的 阵 元 
间距 作为 约束 条 件 ， 选 择 最 优 结 果 ， 确 定 最 终 的 
阵 元 布局 。 

文 [29] 提 出 了 增加 一 个 代价 函数 (Cost Function, 
CF) 的 方法 进一步 限制 优化 结果 ， 即 在 天 空 某 些 
区 域 优化 结果 必须 使 得 代价 函数 最 小 ， 这 样 可 以 
限制 某 些 强 源 对 观测 效果 的 影响 。 但 他 们 的 研究 
并 没有 考虑 稀 布 阵列 的 优势 ， 如 果 能 进一步 降低 
天 线 阵 元 数量 ， 则 可 以 减轻 后 端的 处 理 压力 。 

在 稀 布 方案 方面 ， 文 [30] 分 析 了 规则 稀 布 阵 
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图 4 天 线 阵 列 的 稀 布 + 低 副 办 技术 处 理 后 的 性 能 提升 
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Fig.4 The performance upgrade through sparse and 
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图 5 LOFAR 96 个 LBA 天 线 稀 布 分 布 图 
Fig.5 LOFAR 96 elements sparse configuration 


JJ ( Regular Sparse Array ) 和 非 规 则 稀 布 阵列 (ITregular Sparse Array ) 之 间 的 性 能 对 比 (图 7) ， 通 过 仿真 
得 出 非 规则 稀 布 阵列 具有 更 平 请 的 天 区 扫描 图 像 的 结论 。 原 因 在 于 非 规 则 阵列 的 阵 元 位 置 自 由 度 更 
大 ， 最 优化 搜索 范围 更 大 ， 能 优化 出 更 低 的 最 高 副 办 电 平 和 远 区 栅 辩 电 平 .*1。 但 是 非 规则 稀 布 阵 列 
的 最 优化 搜索 难度 比 规则 稀 布 阵列 大 大 增加 ， 这 是 必须 在 布 阵 之 前 解决 的 问题 。 
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图 6 经 优化 后 LWA-1 天 线 阵 的 布局 情况 。(a) 原始 规则 布局 ; (b) 优 化 后 稀 布 布局 
Fig. 6 Optimized layout of the LWA-1 antenna array. (a) Original regular layout; (b) sparse layout after optimization 
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图 7 规则 稀 布 阵列 和 非 规则 稀 布 阵列 成 图 效果 对 比 


Fig.7 The difference results with regular sparse array 
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一 个 完备 的 宽带 射电 天 文 甚 高 频 天 线 阵 观测 结果 是 多 种 因素 共同 作用 的 结果 ， 主 要 包括 : (1) A 
天 线 阵 元 的 辐射 方向 图 ; (2) 站 址 坐标 ， 包 括 经 纬度 、 海 拔高 度 ; (3) 该 频段 内 天 空 辐射 ， 包 括 多 种 
强 源 的 分 频段 流量 、 干 扰 源 方位 等 ; (4) 各 种 误差 ， 包 括 阵 元 放置 误差 、 幅 相 误差 等 ; (5) 观测 时 间 。 
这 些 因素 的 关系 如 图 8。 
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图 8 决定 天 线 阵 输出 的 多 种 因素 内 在 关系 图 


Fig.8 The relationship of different factors which determine the array output 


如 此 多 的 决定 因素 给 稀 布 阵列 在 射电 天 文 甚 高 频频 段 的 应 用 带 来 了 诸多 技术 挑战 : (1) 稀 布 阵列 
自身 最 优 解 的 搜索 困难 。 稀 布 阵列 的 全 局 最 优 解 是 个 NP 困难 问题 ( Non-deterministie Polynomial-time 
hard, NP-hard) ， 目 前 雷达 、 微 波 领域 仍 未 很 好 地 解决 ， 所 有 优化 搜索 算法 均 陷 入 一 个 局 部 最 优 解 ， 
且 不 同 的 约束 条 件 对 算法 的 要 求 不 同 ， 得 到 最 终 的 布局 方案 也 不 一 样 ， 这 也 是 目前 稀 布 方案 研究 如 火 
如 蔡 的 原因 。 射 电 天 文 的 其 高 频频 段 需 要 考虑 和 雷达 等 领域 不 一 样 的 约束 条 件 ， 目 前 融合 射电 天 文 约 
束 条 件 的 优化 搜索 算法 还 没有 很 好 的 研究 ; 其 高 频 射电 天 文 阵列 的 约束 条 件 究 竞 有 哪些 ， 还 没有 系统 
的 整理 和 研究 。 

(2) 射 电 天 文中 远 区 栅 为 的 影响 。 由 于 射电 天 文 信号 大 多 是 微弱 信号 ,来 自 远 区 的 栅 办 也 干扰 主 
激 观 测 数据 ， 根 据 文 [29] 的 定义 ， 一 个 天 线 阵 列 在 成 图 过 程 中 ， 从 远 副 办 进入 的 噪声 信号 (Far 
Sidelobe Confusion Noise, FSCN) 定义 为 多 个 干扰 源 s gos BE PRO I, 与 天 线 阵 方向 图 的 点 扩散 函数 
( Point Spread Function, PSF)! : 


FSCN = X I, * PSF. (4) 


目前 的 处 理 方案 采用 CLEAN AEO), 3 SERE D Ec eE 69 F DEUS ORG ERE M AEEA 08 RE rg 
号 。 但 是 上 述 方案 的 重要 前 提 是 需要 精确 知道 干扰 源 s 和 点 扩散 函数 之 间 位 置 及 幅度 关系 ， 且 要 求 最 
高 副 办 电 平 要 低 于 主办 30 dB 以 上 ， 才 能 抑制 一 些 强 源 辐射 3 。 

为 此 ,结合 天 线 阵 综合 方向 图 低 副 瓣 技术 , 文 [32] 通 过 仿真 将 泰勒 (Taylor) 加权 和 Doph- 
Chebyshev 综合 加 权 两 种 方式 得 出 结论 ，Doph-Chebysherv 综合 能 抑制 副 淮 电 平 达 30 dB 以 上 ， 这 样 可 
以 有 效 抑 制 大 部 分 强 源 干扰 ,但 是 泰勒 综合 达 不 到 如 此 高 的 抑制 度 。 男 一 方面 ， 文 [31] 只 尝试 了 儿 
种 综合 加 权 算 法 ,希望 后 续 的 研究 能 够 广泛 调研 各 种 综合 加 权 算 法 的 性 能 。 本 文 作者 也 产生 疑惑 ， 在 
阵列 综合 中 ， 泰 勒 综合 及 其 衍生 算法 的 抑制 比 Doph-Chebyshev 在 副 办 抑制 要 好 ， 是 否 是 作者 没有 仿 
真 研究 充分 ? 
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进一步 ， 射 电 天 文 领 域 的 宽带 且 随 机 不 规则 分 布 的 稀 布 阵 ， 更 需要 对 不 同 参数 进行 综合 考量 ， 例 
如 频率 不 同 、 单 元 数 不 同 且 加 权 方 案 不 同时 的 加 权 因 子 ， 这 对 于 稀 布 阵列 的 综合 仿真 来 说 是 巨大 的 挑 
战 。 由 于 其 与 规则 阵列 不 同 ， 要 精确 获知 波束 轮廓 需要 根据 稀 布 布局 规律 进行 仿真 ， 并 且 对 于 宽带 观 
测 的 射电 天 文 阵 列 ， 例 如 跨 3-4 倍 频 程 ， 在 不 同 频段 存在 不 同 的 波束 轮廓 ， 特 别 是 在 多 波段 成 像 中 ， 
要 获知 上 述 多 个 频段 内 的 波束 轮廓 ， 需 要 海量 计算 仿真 。 

(3) 宽带 低频 阵列 的 误差 来 源 多 。 对 于 一 个 单元 数量 巨大 的 天 线 阵 ， 其 引入 的 位 置 、 相 位 和 幅度 
等 误差 对 方向 图 造成 影响 ， 特 别 是 稀 布 阵列 ， 在 位 置 和 幅 相 上 引入 的 误差 更 大 。 现 有 研究 中 的 仿真 基 
本 是 理想 化 的 ， 即 认为 上 述 误 差 不 存 在 ,实际 工作 需要 考量 和 融合 这 些 误 差 。 这 样 的 建 模 分 析 目 前 在 
射电 天 文 领域 还 没有 涉及 ,我 们 需要 考虑 建立 一 个 综合 的 误差 模型 ， 这 样 计算 量 也 是 海量 的 。 

同时 对 多 个 通道 带宽 的 选取 造成 上 下 限 频 率 偏 移 ， 在 仿真 时 一 般 采 用 设置 该 带 内 的 中 心 频 点 为 仿 
真 频 点 ， 获 得 阵 元 加 权 因 子 。 但 是 对 于 该 带宽 内 的 最 高 频率 和 最 低频 率 在 指向 精度 的 偏差 .该 偏差 究 
竞 在 其 他 频 点 上 引起 多 少 误差 ， 目 前 在 均匀 阵 方面 有 一 定 的 研究 ， 但 是 在 宽带 稀 布 阵列 上 还 没有 
相关 的 文献 报道 。 带 宽 的 选择 又 和 前 端的 现场 可 编程 门 阵 列 ( Field Programmable Gate Array, FPGA) 处 
理 能 力 、 系 统 的 传输 能 力 甚 至 数据 存储 能 力 息息相关 。 

(A) 天 空 模 型 和 可 重 构 阵 列 方 向 图 之 间 的 卷 积 耦合 带 来 计算 困难 。 天 空 多 频率 辐射 模型 是 整个 天 
线 阵 特别 是 甚 高 频频 段 天 线 阵 观测 成 败 的 关键 ， 根 据 天 线 阵 亮 温 度 模型 


[| lo, $) PT, C0, $)singdode 
T, = 


ll |F(0, p) | "singdgd 由 


其 中 ， 天 空 背 景 辐射 T (0, 中) 来 源 复杂 ， 既 有 强 射电 源 ， 又 可 能 包含 一 些 干扰 频率 ， 文 [31] 提 到 采 
用 多 种 天 空 背景 模型 精确 计算 Tu,(6, 由 ) ， 但 是 由 于 仿真 平台 的 算 力 ， 还 需要 进一步 深入 开展 这 方面 
的 研究 。 同 时 ， 天 空 模型 是 点 扩散 函数 和 不 同 频段 天 空 背景 模型 的 相互 卷 积 作用 为 进一步 获知 远 副 办 
进入 的 噪声 信号 ， 而 射电 天 空 模型 在 不 同 频率 的 表现 不 同 。 对 于 宽带 天 线 阵 来 说 ,求解 不 同 频率 下 远 
FIERE AIR fme. ， 必 须 考 虑 前 端子 通道 的 划分 ， 以 LOFAR 的 LBA 段 (10~80 MHz) 为 例 ， 经 过 两 
次 信道 化 ， 共 划分 为 512 个 子 通道 ， 每 个 子 通道 带宽 为 195 kHz 或 156 kHz， 这 样 最 好 能 够 在 每 个 子 
通道 上 进行 方向 图 综合 以 及 低 副 办 设计 ， 需 要 巨大 的 算 力 和 通信 能 

文 [29] 根 据 Mort 等 人 提出 的 建议 ， 基 于 更 加 精确 的 射电 天 空 背景 模型 ， 充 分 考虑 银河 系 辐射 、 
多 亮 源 辐 射 的 天 空 分 布 特点 ， 同 时 这 些 辐射 又 与 频率 相关 。 如 果 对 每 个 子 通道 进行 点 扩散 函数 求解 ， 
并 给 出 相应 的 误差 ， 可 以 为 下 一 步 站 站 之 间 的 干涉 测量 中 每 通道 的 精确 指向 、 脏 图 结 化 (CLEAN ) f£ 
序 的 使 用 提供 依据 2 。 多 种 天 空 射电 辐射 分 布 背景 模型 的 使 用 ， 精 分 到 多 个 子 通道 上 的 流量 ， 可 以 
进一步 为 天 线 阵 方向 图 综合 算法 以 及 指向 提供 依据 。 


3 ” 稀 布 阵列 优化 设计 分 析 


3.1 天 线 阵 方向 图 的 一 般 表 达 
一 般 阵 列 方 向 图 综合 问题 可 以 表述 为 


(5) 


E(u) A(T, u,) A(r;, u,) ub ACry, u,) L 
E(u, A(r,, w A(T, us e Á(Ty, u$ L 

ce cms n. cm er] sls (6) 
E(u) Al Tis Uy) A(T, uy) s A(ry, uy) Iy 


其 中 , E=[E(u,), E(u), (uy) ] 是 期 望 的 阵列 方向 图 函数 ; w=[ui, ws,…, uy EARE 
矢量 ; r=[7, 7,,…, rn] 是 阵 元 的 位 置 撩 量 ; 7= [ 1, L, …, W] 是 激励 电流 矢量 。 人 研究 选取 平面 阵 
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列 指向 茶 特 定 方向 (po，b%) 天 线 方向 图 函数 的 普 适 公式 作为 优化 函数 
F(0, e) Š > X0. 8) w, e Ln indes ineo) +d „( singcosp -sin6cosgo) ] : (7) 
进一步 考虑 稀 布 阵列 的 优化 建立 。 首 先 ， 天 线 阵 单元 分 布 -方向 图 综合 模型 的 建立 ， 如 (6) 式 ， 
将 一 个 二 维 平 面 阵 列 在 研究 中 进行 网 格 化 ， 在 阵列 平面 内 切 分 为 mn 个 网 格 ， 即 在 横向 分 为 m íT, 
纵向 分 为 n 列 ， 每 个 网 格 中 至 多 有 一 架 天 线 ， 如 图 9。r, 是 第 mn 个 阵 元 距离 初始 位 阵 元 的 位 置 矢 
量 ， 可 以 分 解 为 x M y 两 个 维度 对 初始 阵 元 的 距离 差 ， 即 [ 4, , 4,,]， 那 么 (7) 式 可 以 改写 为 


N-1 M-1 


F( 0, e) = 2: 2; £. ( 由， O) w el à Ln neos singeong) +dym(singcosp - sinócosgg) ] . (8) 

其 次 ， 阵 列 优化 目标 的 选取 。 根 据 前 述 多 篇 文献 所 述 ， 以 最 高 副 辩 电 平 作为 优化 目标 。 根 据 最 大 

劳 办 电 平 的 定义 ， 取 适应 度 函 数 为 方位 向 方向 图 最 大 旁 兴 电 平 和 俯仰 向 方向 图 最 大 旁 兴 电 平 的 和 |: 

MSLL = max[ Fas( 9, 0,; po, 6,) ] 十 max[ Falo, Oo; $o, 6,) ] 
= a F(o, 0,; po, 0o) | + mad F(o, 0, ; po, 0o) 
pesi max[ F(¢, 05; Po, 0o) ] pes max[ F($, 05; Qo, 8] ' 

其 中 ，max 表示 求 最 大 值 函 数 ;， 5 表示 0 = 0, EZ D 8] Zr T8] PEL DRE DECIR] ,— 575 D 181 2 8] PERS] EIE 

功率 点 为 2p。， 则 


(9) 


Si = {pl Pun pb -po U po tp SP S Pu) ; (10) 
SER e = p, HRES 1681 ELSE CIR] ， 若 俯仰 向 方向 图 的 主办 零 功率 点 为 2， 则 
S, = {0 0an 0x0,-y, UO, * y, x0zx0,.). (11) 
定义 优化 模型 为 
minCMSLL) (12) 


目前 ， 许 多 文献 的 优化 一 般 基 于 遗传 算 
法 ， 但 是 由 于 至 今 没 有 一 个 理论 给 出 稀 布 阵列 
最 高 劳 流 电 平 可 以 达到 的 下 限 ， 为 此 阵列 稀 布 
技术 一 直 以 来 处 于 不 断 改进 和 完善 状态 。 在 遗 
传 算法 基本 流程 的 基础 上 ， 国 内 学 者 提出 了 许 
多 改进 算法 ， 如 针对 传统 遗传 算法 中 “ 轮 盘 
赌 ”方式 选择 下 一 代 染 色 体 导致 算法 收敛 过 
快 ， 优 化 结果 早熟 的 问题 ， 文 [33] 提 出 了 一 种 
改进 的 遗传 算法 (Improved Genetic. Algorithm, 
IGA), 将 当代 最 优 染 色 体 及 具有 较 小 适应 度 
的 前 3 条 染色 体 引 入 下 一 代 种 群 ， 其 余 染 色 体 2 
按照 “ 轮 盘 赌 ”的 方式 选择 。 采 用 这 样 的 方式 ~ Alfy, üy) 
能 够 克服 算法 早熟 ， 并 加 快 算法 的 收敛 速度 。 

文 [34] 提 出 了 一 种 基于 多 种 群 首领 决策 遗 
传 算法 (Multi-Population Leader Dominating Genetic 
Algorithm, MPLDGA), ， 这 种 算法 可 以 避免 传统 的 遗传 算法 在 最 优 搜索 过 程 中 陷 人 局 部 最 优 解 。 但 是 
射电 天 文 天 线 阵 布局 需要 针对 不 同 站 点 、 不 同 观 测 射电 源 性 质 进行 分 类 。 例 如 太阳 这 种 强 射电 源 观 
JW, ， 追 求 较 高 的 时 间 和 频谱 分 辩 率 ， 一 般 太 阳 射 电 爆 发 为 正信 噪 比 信号 ， 甚 至 强 于 部 分 干扰 信号 ， 在 
天 线 阵 布局 优化 方面 对 副 兴 电 平 要 求 相 对 较 低 ;又 如 由 于 宇宙 再 电离 时 期 (Epoch of Reionization ， 
EoR) 和 脉冲 星 等 射电 源 属 于 弱 信 号 ， 甚 至 远 低 于 部 分 射电 源 信号 强度 (最 强 10 000 倍 以 上 ) ， 观 测 中 
需要 长 时 间 积分 ， 在 天 线 阵 布局 优化 方面 需要 更 低 的 副 办 电 平 ， 为 此 文 [32] 建议 ， 副 流 电 平 至 少 要 


图 9 稀 布 阵列 空间 示意 图 


Fig.9 Space Schematic diagram of sparse array 
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低 于 主流 电 平 30 dB, 
3.2 通过 阵 元 加 权 波 束 赋 形 函 数 进一步 实现 低 副 办 

阵列 天 线 的 主要 优点 是 可 以 实现 精确 控制 的 加 权 ， 以 产生 极 低 副 办 的 方向 图 或 非常 接近 选 定形 状 
的 方向 图 。 根 据 天 线 阵 照射 函数 与 副 注 之 间 的 关系 ， 远 区 副 辩 的 衰减 受 口径 场 加 权 函 数 的 导数 控制 。 
我 们 需要 进行 几 项 思考 和 研究 : 首先 ， 拟 尝试 采用 多 种 天 线 阵 综合 算法 ， 对 天 线 阵 方向 图 进行 综合 ， 
同时 对 比 天 线 阵 效率 与 最 高 旁 兴 电 平 抑制 度 之 间 的 关系 。 在 传统 的 研究 中 ， 函 数 及 其 导数 中 具有 不 连 
续 点 的 均匀 加 权时 ， 远 区 副 洲 随 (Zu/A) 一 变化 。 如 在 采用 三 角 函 数 加 权时 ， 远 区 副 办 随 (Ze 人 A) 一 变 
化 ， 而 采用 余弦 平方 加 权时 ， 远 区 副 办 随 ( Ze/A) 一 变化 。 

对 于 复杂 问题 (例如 多 目标 优化 问题 和 同时 优化 多 种 阵列 参数 问题 等 ) 难 以 找到 确定 的 综合 公式 。 
随 着 计算 机 技术 的 发 展 和 阵列 天 线 方 向 图 性 能 要 求 的 不 断 提高 ， 如 今 多 采用 优化 方法 解决 综合 问题 。 
随 着 低 副 办 技术 的 深入 研究 ， 采 用 不 同 的 解析 多 项 式 结 构 实 现 加 权 函 数 ， 根 据 天 线 主办 宽度 、 最 高 责 
办 电 平 等 约束 条 件 ， 求 解 每 个 阵 元 每 个 频 点 的 加 权 人 参数 。 

近年 来 各 种 优化 算法 纷纷 涌现 ， 主 要 有 Chebyshev 综合 、 泰 勒 综合 、 储 里 叶 变 换 法 、Woodward 
综合 技术 、 交 错 投影 、 功 率 方向 图 综合 法 等 一 系列 方法 ， 非 常 适 用 于 波束 赋 形 。 和 常用 的 有 Chebyshev 
综合 和 泰勒 综合 及 其 改进 型 综合 法 。 射 电 天 文 甚 高 频 天 线 阵 研究 将 比 对 多 种 算法 的 赋 形 效果 ， 特 别 是 
针对 射电 天 文中 的 特殊 要 求 ， 在 副 泊 抑制 度 、 波 束 稳定 性 和 灵敏 度 变 化 之 间 求 取 平 衡 。 

其 次 ， 对 于 不 同 的 指向 角度 (9, p), RER A 时 的 不 同 波束 合成 方案 的 效果 情况 ， 该 效果 的 
评价 应 为 射电 天 文 需 求 为 导向 ， 考 察 如 下 几 点 (1) 考察 整体 天 线 阵 的 最 高 旁 闪 电 平 在 指 问 变化 后 的 
了 恶化， 特别 是 在 某 些 俯仰 角 以 下 时 ， 最 高 劳 兴 电 平 难以 用 加 权 综 合 的 方法 进行 抑制 ， 依 此 确定 是 否 为 
该 天 线 阵 的 最 佳 可 观测 天 区 ; 根据 前 述 不 同 射电 源 观 测 需 求 ， 优 化 选择 不 同 的 加 权 函 数 ， 特 别 是 不 同 
观测 目标 对 不 同 副 辩 电 平和 波束 赋 形 带 来 的 增益 变化 等 需求 程度 进行 比 对 分 析 。(2) 融 合 已 有 的 其 高 
频频 段 这 天 数据 ， 研 究 天 空 背 景 辐射 分 布 的 融合 介入 绽 合 布局 效果 分 析 ， 通 过 亮 源 分 布 ， 分 析 天 线 阵 
综合 后 的 点 扩散 函数 与 天 空 背景 的 卷 积 效 果 ， 分 析 此 时 不 同 频段 的 站 波束 稳定 性 。 

3.3 多 误差 模型 来 源 分 析 

天 线 阵 方向 图 的 好 坏 取决 于 多 种 误差 ， 这 些 误差 在 目前 的 射电 天 文 天 线 阵 设计 中 人 鲜 有 考虑 到 实际 
的 方向 图 综合 ， 主 要 有 幅度 误差 、 相 位 误差 和 加 权 频 率 偏差 等 。 

3.3.1 阵列 幅度 误差 

阵列 的 幅度 误差 主要 由 阵 元 天 线 及 接收 机 增益 响应 不 一 致 导致 的 ， 大 多 数 由 环境 温度 、 接 收 机 幅 
度 稳定 性 引起 ， 一 般 认 为 幅度 误差 与 信号 的 来 波 方向 无 关 。 幅 度 误 差 是 天 线 阵 元 因子 在 幅 值 上 引入 的 
一 个 乘 性 误差 .通过 两 两 天 线 之 间 的 幅度 差 测 量 进行 确定 。 我 们 拟 研究 测量 值 的 误差 范围 ， 并 建立 幅 
度 误差 模型 。 

如 (7) 式 ,假设 一 个 m xn 元 阵列 ， 理 想 条 件 下 阵 元 的 天 线 方向 图 f,($, 9) ， 当 阵列 存在 幅度 误 
差 时 ， 实 际 的 导向 矢量 为 


Fab, 0 =T, f. (o, 0), (13) 

其 中 , T lp, py ns pu ] 为 阵列 的 幅度 误差 矩阵 ，p RREFERA ij 列 的 阵 元 的 幅 

度 误 差 (不 一 定 是 规则 分 布 ) 。 在 研究 中 将 接收 机 的 误差 折算 到 天 线 方向 图 中 ， 则 方向 图 模型 可 以 改 
写 为 

F(80, Q) - * > ouo. 0)w " el à. Ln sinBcose -sinfgcosg) * d, ( sin&cosp -sincosQo) ] . (14) 


3.3.2 阵列 相位 误差 

阵列 的 相位 误差 主要 由 阵 元 天 线 及 接收 机 相位 啊 应 不 一 致 导 致 的 ， 大 多 数 由 环境 温度 和 接收 机 相 
位 稳定 性 引起 ， 一 般 认 为 相位 误差 与 信号 的 来 波 方向 无 关 。 相 位 误差 是 天 线 阵 元 因子 在 相位 上 引入 的 
一 个 加 性 误差 ， 需 要 通过 两 两 天 线 直 线 的 相位 差 测 量 进 行 确 定 。 我 们 拟 研 究 测量 值 的 误差 范围 ， 并 建 


430 X x 研究 与 d 20 卷 


立 相 位 误差 模型 。 如 (7) 式 ,假设 一 个 mxn 元 阵列 ， 阵 元 的 位 置 矢量 为 r=[7 ra, s Tm], EA 
存在 相位 误差 时 ， 实 际 的 导向 矢量 为 


一 一 


f(D, OSL, fo, 0), (15) 


KB, falb, 9) 为 理想 条 件 下 的 阵列 导向 矢量 ; Dlh, eft, s, et ] 为 阵列 的 相位 误差 矩阵 ; 
中 ;表示 处 于 二 维 平面 中 第 i 行 j 列 阵 元 的 相位 误差 (不 一 定 是 规则 分 布 ): 


N-1 M-1 


F( 0 op ) 2 2. Mf ( 中 0) p Six [ d, ( sincosg -sin8gcose) *d,,( sinfcosp -sincoseo) ] tonm? ( 16 ) 
, mn , mn . 
n=0 m=0 


3.3.3 加 权 频 率 偏差 

加 权 频 率 偏差 是 重要 的 仪器 误差 效应 ， 在 超 宽带 天 线 阵 中 ， 数 字 阵 列 可 以 通过 对 阵 元 接收 的 不 同 
频率 的 信号 加 相位 权 以 形成 特定 的 指向 波束 。 但 在 实际 应 用 中 ， 系 统 通常 无 法 对 每 个 频率 的 信号 添加 
对 应 的 权 值 ， 而 是 采取 一 定 带宽 的 信号 按照 中 心 频率 加 权 ， 这 就 产生 非 中 心 频率 处 的 加 权 偏 差 ， 需 要 
建立 误差 模型 进行 分 析 。 

针对 前 端 预 处 理 现 场 可 编程 门 阵 列 进行 的 信道 化 处 理 形成 的 中 心 频率 上 万， 带宽 为 BW, ERU 
为 +， 下限 频率 为 及 ~ 。 将 上 述 两 个 频 点 代入 方向 图 和 天 线 阵 进行 综合 计算 ， 重 新 获得 天 线 


方向 图 ， 并 得 到 副 辩 电 平 的 变化 情况 ， 即 AMSLL 以 及 最 高 旁 兴 电 平 处 的 增益 变化 情况 AGysj。AGys 
的 恶化 程度 是 BW 的 函数 ， 在 研究 中 我 们 拟 找 到 在 不 同 布局 、 综 合 情 况 下 的 函数 关系 

AC yj; - F(BW), (17) 
并 在 最 终 的 方向 图 计算 中 加 以 误差 说 明 。 模 型 建立 如 下 : 假设 某 阵列 一 个 数字 通道 的 带宽 为 BW, W 
系统 加 权 后 的 最 大 频率 偏差 为 BW/2, WEN fua. c 


fa -4, (18) 
H 
则 
F ENa. (19) 
L 
引入 加 权 频 率 偏差 后 的 阵列 方向 图 函数 可 以 表示 为 
N-1 M-1 (s } 
F(0 ; - ， 0 Mo OsaND [ d, ( sincosg —sin6gcose ) +d ( sinfcosp -sin&coseo) | ] 20 
(0, o) 2, f. (0 ) wine (20) 


其 中 ,NN 为 整个 观测 频段 上 限 频率 所 与 该 信道 中 心 频率 的 比值 。 
3.3.4” 阵 元 位 置 误差 

天 线 阵 列 由 于 安装 工艺 等 影响 不 可 避免 地 存在 阵 元 位 置 误差 ， 阵列 综合 给 出 的 各 阵 元 加 权 参 数 相 
位 是 基于 阵 元 的 理想 位 置 计算 确定 的 ， 所 以 阵 元 的 位 置 误 差 导致 阵列 方向 图 的 变化 ， 并 有 可 能 导致 副 
办 电 平 比 仿真 结果 有 所 提高 。 

在 研究 中 我 们 拟 建立 一 个 阵 元 位 置 误差 范围 ， 在 方向 图 计算 时 考虑 该 误差 范围 的 上 下 限 恶 化 程 
度 ， 并 在 最 终 的 方向 图 计算 数据 中 加 以 误差 说 明 ， 模 型 建立 如 下 : 如 (7) 式 ,假设 一 个 m xn 元 阵列 ， 
阵 元 的 位 置 矢 量 为 rz=[m , ros. 7,,] ， 当 阵列 存在 位 置 误差 时 ， 实 际 位 置 矢量 为 

r'-[r + Ar, r + Ar, 0, Ta + Anl, (21) 


其 中 ，Ar, 可 以 分 解 为 Ad, 和 Ad,， 则 引入 阵 元 位 置 误差 后 的 阵列 方向 图 函数 可 以 表示 为 


N-1 M-1 2n 


F( 0 o ) = > ` f. ( 中 8) w el n [ (d, * Ad, ) ( sin&cosp -singcose) * ( dj, * Ad,,) ( sincosp —sinbcosp0) ] ( 22 ) 
? mn E mn š 
n=0 m=0 
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3.3.5 综合 偏差 模型 的 建立 

上 述 4 个 误差 模型 均 只 考虑 单 种 误差 的 情况 ， 在 实际 的 天 线 阵列 中 各 种 误差 共同 存在 ， 最 终 研 究 
将 建立 一 个 综合 方向 图 误差 引入 模型 。 

上 述 各 种 误差 模型 均 是 随机 误差 的 复 指数 函数 ， 对 于 均匀 阵列 而 言 ， 可 以 用 数学 解析 方法 近似 获 
得 最 高 劳 办 电 平 ， 但 对 于 阵 元 非 均匀 排列 的 稀 布 阵 列 而 言 ， 很 难 用 数学 解析 的 方法 找到 副 准 电 平 出 现 
的 位 置 并 求解 最 高 旁 欠 电 平 。 

通过 计算 机 模拟 手段 ， 可 以 在 已 定义 好 的 稀 布 阵列 因子 仿真 计算 程序 中 融入 多 个 误差 模型 ， 通 过 
测量 + 预 佑 等 方法 ， 设 定 误差 的 上 下 限 ， 分 别 计 算 可 能 的 恶化 程度 。 


4 信和 号 处 理 体系 构建 


4.1 系统 构架 

目前 我 们 已 经 得 到 第 1 步 的 经 费 支 持 ， 建 设 
中 ( 国 ) 马 (来 西亚 ) 两 地 的 其 高 频频 段 天 线 阵 ， 
前 期 开展 了 天 线 单元 、 采 集 系 统 和 数据 预 处 理 方 
面 的 研究 ， 在 已 有 的 16~32 阵 元 其 高 频 天 线 阵 
上 开展 侧重 于 稀 布 阵列 布局 、 综 合算 法 以 及 抗 干 
扰 等 方面 的 研究 ， 并 以 该 频段 重要 的 太阳 射电 观 
测验 证 为 第 1 步 实验 观 测 步骤， 可 以 将 前 期 各 方 
面 的 研究 成 果 应 用 于 该 实验 阵列 。 

借鉴 AAVS( Aperture Array Verification System) 
阵列 的 经 验 ， 前 期 实验 阵列 采用 较 少 阵 元 在 天 线 
阵 元 、 接 收 机 、 数 字 采 集 、 波 束 合成 信号 处 理 等 


方面 进行 实验 性 研究 ， 如 图 10 “”， 为 后 续 大 图 10 AVVSO 用 于 验证 稀 布 阵列 的 16 单元 阵列 

规模 阵列 部 署 扫 清 技术 障碍 。 Fig. 10 AVVSO used to verify a 16-element array 
首先 ， 根据 台 站 的 综合 条 件 ， 如 天 线 阵列 布 in a sparse array 

局 尺寸 的 显示 、 阵 元 间距 (最 大 、 最 小 间距 ) 以 及 


无 线 电 环境 等 ， 采 用 不 同 算法 遵 选 最 优 的 天 线 阵 布 局 方案 ; 第 2 步 ， 在 已 确定 阵列 布局 的 条 件 下 ， 通 
过 各 种 天 线 加 权 算 法 对 不 同 频率 、 不 同 扫 描 角 度 进行 波束 赋 形 计算 ， 并 得 到 加 权 因 子 ; 第 3 步 ， 通 过 
网 络 向 前 端 负责 信号 预 处 理 的 现场 可 编程 门 阵列 进行 加 权 因 子 下 发 ; 第 4 步 , 代入 前 述 各 种 误差 模型 
计算 得 到 不 同 频率 时 的 天 线 阵 方向 图 、 灵 敏 度 等 性 能 参数 ， 最 后 ， 通 过 实时 多 频 天 空 背景 模型 与 上 述 
天 线 阵 卷 积 ， 通 过 快 视 等 手段 评估 当前 观测 效果 。 综 上 ， 该 阵列 的 系统 处 理 流程 如 图 11。 
4.2 ”硬件 平台 的 搭建 

系统 硬件 如 图 12。 硬 件 采用 图 形 处 理 器 或 云 计 算 服 务 器 一 运算 卡 一 高 速 网 络 通信 一 SOPC 协议 
分 发 一 控制 高 速 模 数 转换 器 采集 架构 。 首 先 ， 采 用 图 形 处 理 器 或 者 云 计 算 服 务 器 ， 对 天 线 阵 全 频段 
内 的 多 个 信道 进行 方向 图 综合 计算 、 误 差分 析 、 数 据 传 输 与 存储 等 工作 ， 这 是 典型 的 并 发 计算 。 为 
此 ， 图 形 处 理 器 内 的 计算 需要 考虑 的 问题 较 多 ， 如 综合 误差 模型 、 多 频 点 方向 图 的 计算 等 。 对 于 超 宽 
带 阵 列 ， 我 们 需要 对 每 个 子 通道 进行 单独 分 析 ， 这 样 的 计算 量 是 巨大 的 。 利 用 图 形 处 理 器 实现 多 线程 
的 高 速 编程 ， 且 要 和 前 端的 信道 化 数量 匹配 ， 达 到 最 佳 的 计算 效果 ， 人 快速 计算 分 析 每 个 子 通道 方向 
图 ， 以 及 在 误差 加 扰 时 的 恶化 程度 ， 进 一 步 计 算 该 通道 (中心 频 点 ) 的 加 权 参 数 。 全 频段 内 的 多 个 信 
道 进行 天 线 阵 方向 图 综合 计算 、 误 差分 析 、 数 据 传输 与 存储 等 工作 。 然 后 ， 前 端 采集 和 信和 号 处 理 采 用 
多 通道 ADC + SOPC 的 现场 可 编程 门 阵列 架构 。 通 过 有 效 的 SOPC 配置 ， 实 现 图 形 处 理 器 计算 平台 到 
前 端 现场 可 编程 门 阵列 预 处 理子 系统 的 快速 通信 、 配 置 ， 将 不 同 通道 的 阵列 因子 参数 利用 SOPC 进行 
高 速配 发 。 
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图 11 系统 处 理 方案 流程 图 
Fig. 11 Flowchart of system processing pipeline 
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图 12 系统 硬件 图 
Fig. 12 System hardware diagram 


4.3 基于 图 形 处 理 器 平台 的 多 通道 天 线 阵 列 因 子 及 角度 -灵敏 度 计算 


FPGA 
board-1 


FPGA 
board-2 


FPGA 
board-N 


图 形 处 理 器 岁 形 化 编程 平台 具有 单 指令 多 线程 (Single Instruction Multiple Threads, SIMT) 运行 模 
式 ， 因 而 非常 适用 于 数据 并 行 计算 。 首 先 在 图 形 处 理 器 内 进行 分 析 计 算 ， 然后 通过 GPU-FPGA 的 数 
据 传输 分 发 。 阵 列 方向 图 在 不 同 频率 下 计算 甚 高 频频 段 宽 带 阵 列 ， 如 SKA-low 预计 从 70 MHz 一 直到 
350 MHz， 相 对 带宽 达到 5 : 1， 在 高 速 采集 后 ， 现 场 可 编程 门 阵列 内 部 进行 信道 化 处 理 ，LOFAR 初步 
将 观测 带宽 分 为 512 个 子 带 ， 每 个 子 带 195 kHz( 100 MHz 带宽 ) 或 156 kHz( 80 MHz 带宽 ) 。 根 据 前 面 


天 线 方向 图 综合 (7) 式 ， 微 波 学 中 定义 一 个 特定 波长 A; 电磁 波 空间 波 数 为 
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lii 


pay 


k; B A , (23) 
针对 一 个 频 点 将 (7) 式 改写 为 
N-1 M-1 
F.(0, Q) E ba NGC B) w, ei Leo infos ingoeovg) * d, ( sinfcosg -sinfcosgp) ] . (24) 
n=0 m=0 


进一步 考虑 男 外 两 个 参数 ， 即 每 个 阵 元 在 不 同 频 率 时 的 响应 f, (b, 9) ， 以 及 不 同 频率 时 的 加 权 因 子 ， 
则 应 改写 为 


foa, 0) fuo, 8), (25) 
表示 位 于 二 维 网 格 点 中 ， 第 mn 个 网 格 中 的 天 线 单 元 在 波长 A; 时 方向 图 ; 
Ven? Mi (26) 


表示 位 于 二 维 网 格 点 中 ， 第 mn 个 网 格 中 的 天 线 单元 在 波长 A; 时 加 权 因 子 ， 则 (7) 式 修正 为 


F( 0 Q) = > by { 中 0) w girl da sinteosp-sinfgosp) +dm(singcosp -sin6coseo) ] ( 27) 
i 7 mn,i , mn,i i 


根据 天 线 阵 方向 图 等 于 阵列 因子 和 天 线 因子 的 乘积 ， 以 FCO, o) 加 单元 天 线 的 增益 (设计 保证 ) , 
就 是 目前 在 (9，9p) 方 向 时 的 天 线 阵 增益 G(0, p)。 根 据 公式 求 得 天 线 阵 此 时 的 有 效 面积 > 


4TC(O， 
AG, 网 = 0) 


根据 (7) 式 可 以 计算 得 到 当前 的 灵敏 度 。 充 分 利用 图 形 处 理 器 的 多 线程 并 行 同步 计算 能 力 ， 将 (27) 式 
中 不 同 频 率 的 综合 加 权 因 子 、 方 向 图 、 有 效 面积 以 及 灵敏 度 等 参数 进行 计算 ,流程 如 图 13。 
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图 13 图 形 处 理 器 内 的 计算 流程 图 
Fig. 13 Computing pipeline in CPU 
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本 文 需要 研究 的 主要 问题 如 下 : (1) 如 何在 图 形 处 理 器 内 多 个 线程 并 行 计算 不 同 频 率 时 的 天 线 阵 
加 权 参 数 ， 进 一 步 在 获得 加 权 参 数 后 计算 得 到 此 时 不 同 频率 的 天 线 阵 方 向 图 、 有 效 面积 和 灵敏 度 。 在 
综合 加 权 参 数 的 计算 中 ， 需 要 和 前 端的 天 线 数量 N,、 输 出 极 化 数量 N, 以 及 每 个 天 线 阵 元 输出 的 信道 
化 数量 N, 相 匹 配 ， 这 就 需要 对 图 形 处 理 器 内 可 以 实现 的 线程 计算 能 力 进 行 分 析 。 方 向 图 计算 需要 和 
前 端的 信道 化 数量 N, 相 匹 配 ， 即 图 形 处 理 器 的 线程 处 理 能 力 要 并 行 同时 完成 上 述 IN, 个 信道 频段 的 方 
向 图 。 

(2) 多 线程 计算 的 同步 性 问题 。 在 多 线程 计算 时 还 必须 考虑 它们 之 间 的 同步 性 问题 ， 避 免 由 于 时 
间 不 同步 带 来 的 阵 元 加 权 偏 差 ， 为 此 需要 多 个 线程 内 程序 的 运行 时 间 稳 定 且 保持 一 致 。 

(3) 多 外 部 数据 库 协 调调 用 问题 。 一 些 常 用 数据 来 自 多 个 数据 库 : 志 单 天 线 阵 元 不 同 频率 方向 图 
数据 库 ， 即 f,, (5, 9) ， 在 不 同 频率 时 的 天 线 单元 方向 图 ， 这 可 以 由 HFSS 和 CST 等 软件 获得 一 致 性 
较 好 的 数据 ; @ 综 合 得 到 的 加 权 参 数 数据 库 ， 即 w,, ;， 在 不 同 频率 时 每 个 阵 元 的 加 权 参 数 ， 这 是 由 
上 述 第 2 个 研究 方面 通过 不 同 的 综合 方法 得 到 的 ; 全 射电 源 分 布 数据 库 ， 采 用 上 述 巡 天 数据 ， 配 合 地 
理 经 纬度 、 时 间 等 信息 ， 可 以 快速 计算 得 到 在 天 线 阵 视 场 内 的 位 置 ; 由 多 误差 参数 数据 库 ， 这 是 根据 
天 线 阵 实际 测量 得 到 数据 的 误差 范围 确定 的 。 在 图 形 处 理 器 计算 中 ， 根 据 划分 的 通 带 数 启动 多 个 线 
程 ， 每 个 线程 分 别 调用 上 述 数据 库 中 相应 频率 的 数据 。 

4.4 可 行 性 分 析 
4.4.1 图 形 处 理 器 计算 平台 的 可 行 性 

图 形 处 理 器 计算 平台 具有 强大 的 计算 能 力 ， 目 前 已 经 广泛 应 用 于 天 文 领域 ， 在 脉冲 星 搜 寻 、 大 规 
模 图 像 处 理 等 领域 均 有 应 用 。 因 为 本 项 研究 拟 实现 的 平台 中 ， 图 形 处 理 器 仅 承担 对 阵列 因子 的 计算 以 
及 传输 ， 方向 图 的 生成 和 存储 。 根 据 波 束 宽度 、 地 面 站 与 观测 源 之 间 的 相对 运动 速度 ， 在 主 波束 较 宽 
的 情况 下 ， 观 测 终 端 并 不 需要 非常 短 的 时 间 更 新 一 次 阵列 因子 。 

我 们 拟 采 用 的 图 形 处 理 器 配置 gs4820 服务 器 ,支持 8 路 GPU 运算 ， 具有 8 个 PCI-E * 16 fdv, 
支持 挂 接 24 个 硬盘 ， 最 高 可 扩展 至 3TB 内 存 ， 集 成 了 双 千 兆 网 卡 。 运 算 卡 采用 NVIDIA TESLA V100 
32G 运算 显卡 ， 具 有 5 120 个 运算 核心 (CUDA Cores) ， 在 双 精 度 下 具有 7 Teraflops ， 单 精度 下 具有 14 
Teraflops 的 加 速 性 能 ， 内 存 容量 32 GB ， 显 存 带 宽 为 900 GB/s。 

根据 现 有 的 天 线 阵 体 量 ， 按 照 上 千 架 天 线 和 每 个 天 线 产生 的 上 千 通道 的 数据 流 进 行 分 析 ， 假 设 天 
线 阵 元 数量 为 1 000 架 ， 双 极 化 输出 ， 预 处 理 生 成 的 子 通道 数 为 1024， 每 个 通道 加 权 参 数量 化 位 数 为 
4 Byte。 天 线 阵列 跟踪 精度 决定 主 波束 指向 每 秒 更 新 1 次 ， 采 用 多 波束 观测 时 ， 一 共产 生 8 个 主 波束 ， 
则 每 秒 产生 的 计算 数据 量 为 

D = 1 000 x 2 x 1024 x 4 x 8 = 65.536 MB/s = 560 Mbps . (29) 

首先 ， 服 务 器 内 存 最 大 可 以 集成 3 TB ， 同 时 根据 PCI-E * 16 的 通信 能 力 ， 服 务 器 能 够 完成 传输 任 
务 。 另 一 方面 ， 计 算 卡 需要 开辟 至 少 2 048 个 线程 ， 同 时 每 秒 计 算 量 大 约 为 70MFLOPS ， 这 是 我 们 选 
择 的 图 形 处 理 器 计算 卡 能 够 承受 的 。 

4.4.2 SOPC 内 实现 通信 和 数据 分 发 的 可 行 性 

高 性 能 现场 可 编程 门 阵列 内 的 SOPC 片 可 以 编程 操作 系统 ， 具 有 编程 灵活 方便 等 优势 。C 语言 编 
程 同 时 可 以 封装 自己 的 IP 内 核 , 采用 这 种 系统 很 容易 对 通信 数据 进行 解码 和 分 发 。 

SOPC 采用 基于 ZYNQ XC7Z045 FFG 900 的 高 性 能 计算 板 卡 ， 具 有 高 性 能 ZYNQ-7000FPGA ; 
XC7Z045， 集 成 了 高 速 10-100-1 000 Mbps Ethernet (RGMII) 高 速 网 络 资源 ， 用 于 与 图 形 处 理 器 的 高 
速 通信 ;同时 具有 千 万 门 级 的 逻辑 资源 ， 实 现 对 高 速 ADC 采集 的 控制 、 信 道 化 和 加 权 。 

在 数据 传输 与 处 理 中 ， 主 频 为 2. 5DMIPS(2. 5x 100 万 条 指令 每 秒 ) 的 ZYNQ-7000 系列 FPGA 完全 
能 够 承担 560 Mbps 这 样 的 计算 工作 (ZYNOQ 千 兆 网 口 速率 测试 结果 为 744 ~ 839 Mbps, htips://blog. 
csdn.net/fzktongyong/ article/ details/ 83855974 ) 。 


5 期 


5 结 


n 
T 


亮 等 : 射电 天 文 30~300 MHz 频段 稀 布 阵列 关键 技术 分 析 与 研究 435 


lix; 
T 


it 


其 高 频频 段 射 电 天 文 要 实现 良好 的 观测 性 能 ， 取 决 于 多 种 技术 的 相互 配合 以 及 台 站 条 件 等 方面 的 
综合 考量 。 本 文通 过 对 其 高 频 射 电 天 文 天 线 阵 的 技术 梳理 ， 提 出 了 针对 实验 阵列 的 技术 构架 和 研究 思 
路 。 接 下 来 该 成 果 将 应 用 于 第 1 套 全 数字 化 其 高 频 射电 天 文 阵 列 ， 并 在 试验 中 进一步 调整 和 细 化 系统 
结构 ， 以 期 达到 优良 的 观测 性 能 。 
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Abstract; VHF ( Very High Frequency) band is one of important observation bands in Radio astronomy 
field. In this band, we usually use the array to form phased array telescope. In another hand, sparse array has 
more space resolution and lower side lobe level comparing with the regular array. Furthermore, based on array 
integrated optimization algorithms, we can effectively form the main beam and reject the Maximum side lobe 
levels and far side lobe levels. In this paper, at first, we review the development and research status of VHF 
sparse array in radio astronomy filed and analyze difficulties. The second, we promote a pipeline for sparse 
array building. The first step is constructing the optical elements configuration based on some optimization 
algorithms. And then, we promote a signal processing structure based on high performance computing severs 
and FPGA (Field Programmable Gate Array) SOPC (System On Programmable Chip) platform. In this platform, 
taking advantage of strong computing capability, in the GPU embedded in severs can compute the weighting 
parameters of different elements and different frequency channels. And then, transmitting these parameters 
through the high-speed bus and delivery them to each element in FPGA SOPC. We also analyze the data rate 
with multi-beams that can transmit the weighting parameters in real-time. This work can provide a technical 
base for future large-scale VHF array building. 
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